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In a Mossbauer experiment, if a spa- 
tially-extended absorbing sample is ro- 
tated relative to a moving y-ray source, 
lines of constant y-ray Doppler shift are 
generated through the absorber parallel 
to the motion of the source. As a result, 
resonant absorption takes place along a 
series of parallel lines cutting through 
the absorber, where a particular line is 
determined by the velocity of the 
source. The result is a series of measure- 
ments of line integrals of the absorption 
coefficient through the absorber. An im- 
age or spatial map of the absorption co- 
efficient distribution may then be 
reconstructed using tomographic image- 
reconstruction algorithms. Moreover, 



when measurements are recorded both 
as a function of the source velocity and 
the absorber rotational velocity, spectral 
information may also be recovered as a 
function of position. Spatial resolution is 
proportional to the rate of rotation of 
the absorber, but is ultimately signal-to- 
noise limited. 
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1. Introduction 



The Mossbauer effect is an important spectro- 
scopic tool in materials science [1,2]'. In a typical 
Mossbauer experiment, y-rays emitted by a ra- 
dioactive source are allowed to impinge upon a 
sample of absorbing material; the transmitted or 
resonantly-scattered y-rays are then subsequently 
detected (fig. 1). Thus, in a transmission experi- 
ment, a detector counts the number of y-rays pass- 
ing through the absorber, and in a scattering 
experiment the detector counts the y-rays re-emit- 
ted after resonant absorption. The Mossbauer ef- 
fect, or nuclear y-ray fluorescence, is made 
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possible by the recoilless emission and absorption 
of the y-rays in nuclei embedded in a solid lattice. 
As a result, the resonant linewidths, relative to the 
y-ray energy, are characteristically very narrow. 
Moreover, the resonant absorption can be "tuned" 
merely by moving the source at a small velocity 
relative to the absorber, which imparts some ki- 
netic energy, or Doppler shift, to the y-ray. 
Changes in source-absorber speeds as small as a 
few millimeters per second are often sufficient to 
destroy the resonance. In the conventional Moss- 
bauer experiment, an absorption spectrum of the 
material under study is generated by moving the 
source relative to the absorber and counting the 
transmitted y-rays (or resonantly-scattered y-rays 
in a scattering experiment) as a function of the rela- 
tive source-absorber velocity. 
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Figure 1-Conventional Mossbauer experimental arrangement 
with stationary absorber and moving source. 

In this paper, we propose the idea of Mossbauer 
imaging. An image is by definition some quantity 
or parameter of interest displayed as a function of 
position, i.e., a picture of the parameter as it is dis- 
tributed in space. In Mossbauer imaging, the sim- 
plest example of such a parameter is the Mossbauer 
absorption coefficient, and the aim would be to re- 
construct and display its spatial distribution. A con- 
ventional Mossbauer experiment measures only a 
bulk average of the resonant absorption coefficient 
over the absorbing specimen. Thus, spatial inhomo- 
geneities within an extended absorber, due either to 
variations in the absorption cross section of a pure 
material or to the mixture of different nuclei within 
a composite material, are inevitably averaged out 
in the bulk measurement process. A Mossbauer 
imaging experiment, however, would permit the 
reconstruction of a two-dimensional map of the 
Mossbauer absorption coefficient, that is, a two-di- 
mensional picture of the strength of the resonant 
absorption as a function of position. In another, 
somewhat more complex, version of an imaging ex- 
periment, spectroscopic information as a function 
of position should also be recoverable, rather than 
merely differences in the absorption coefficient as a 
function of position. As a consequence, in the latter 
version, true Mossbauer spectroscopy can be per- 
formed in an imaging mode, so that heterogeneous 
or composite samples may be investigated. 

The idea of Mossbauer imaging was inspired by 
the success of nuclear-magnetic-resonance (NMR) 
imaging, since NMR and the Mossbauer effect 
share some fundamental characteristics, both being 
nuclear resonance phenomena. NMR imaging has 
recently found notable success in diagnostic 
medicine [3]. While there are no foreseeable appli- 
cations of Mossbauer imaging in medicine, applica- 
tions to materials science are thought to exist. The 
ability to perform Mossbauer spectroscopy in an 
imaging mode, i.e., to do spectroscopy as a func- 
tion of position within a sample, rather than in 
bulk, should prove to be of value in the analysis of 
heterogeneous materials. Although NMR and 
Mossbauer imaging may have different applica- 



tions, it is instructive to compare the two tech- 
niques. In NMR imaging, resonant-frequency 
information is translated into spatial information by 
imposing a magnetic field gradient on a system of 
precessing nuclear spins. In Mossbauer imaging, 
we will see that by rotating the y-ray absorber, a 
velocity gradient is imposed along the radial direc- 
tion of the absorber and the resonance associated 
with a particular Doppler shift (analogous to a par- 
ticular NMR precessional frequency) takes place 
along a family of lines perpendicular to the veloc- 
ity gradient (analogous to the NMR magnetic field 
gradient). This process gives rise to a set of mea- 
surements of the line integrals of the absorption co- 
efficient, from which an image can be recovered 
using well-known tomographic image reconstruc- 
tion algorithms. These algorithms are mathemati- 
cally similar to those employed in both x-ray and 
NMR tomography [4,5]. A more detailed descrip- 
tion of the process by which the line-integral mea- 
surements arise is given in section 3. 

The basis of the imaging approach can be most 
easily described using a classical interpretation of 
Mossbauer absorption, which will suffice for our 
purposes. As noted earlier, in this description the 
relative source-absorber velocity imparts a 
Doppler shift to the y-ray. For purposes of illustra- 
tion, consider the simplest case of a material with a 
single resonance line (neglecting isomer shifts); 
then resonance will take place at zero Doppler 
shift, or zero relative velocity between the source 
and absorber. By rotating an extended absorbing 
object relative to the source, lines of constant 
Doppler shift (or equivalently, of constant y-ray 
energy) are generated across the absorber. As a re- 
sult, resonant absorption takes place along one line 
at a time, giving rise to one line-integral measure- 
ment of the absorption coefficient. The location of 
the line is determined by three parameters: the ab- 
sorber rotational velocity, the source velocity and 
the instantaneous position of the absorber during 
the measurement. From a complete set of such line- 
integral measurements, a spatial map of the y-ray 
absorption can, in principle, be tomographically re- 
constructed. We shall see that the spatial resolution 
within the reconstructed image is proportional to 
the ratio of the natural linewidth of the resonance 
(in units of velocity) to the rate of rotation of the 
absorber, but is ultimately signal-to-noise limited. 
This means that, in principle, arbitrarily high reso- 
lution can be achieved, but at the expense of 
rapidly increasing detector integration times. In 
particular, the integration time will increase in di- 
rect proportion to the number of resolution ele- 
ments (pixels) in the desired image. Thus, in a 
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two-dimensional image, doubling the linear resolu- 
tion will square the number of pixels, which in turn 
will square the required integration time. As a con- 
sequence, very high resolution will require exceed- 
ingly long integration times unless sources of high 
intensity become more readily available. 

At this stage one might observe that the equiva- 
lent line integral measurements could be generated 
merely by directing a collimated beam of y-rays at 
the absorber over a wide range of incident angles. 
This approach is analogous to conventional x-ray 
computerized tomography in a transmission experi- 
ment or to single photon emission tomography in a 
scattering experiment. The advantage of the pro- 
posed method is that an extended source can be 
used rather than a collimated source, which implies 
a substantially greater incident y-ray flux and a 
corresponding enhancement in the signal-to-noise 
ratio. 

In the following three sections, we consider the 
idealized case of the y-rays incident parallel to the 
x-axis. This can be achieved by moving the source 
sufficiently far away so that the source subtends a 
small angle at the absorber. The assumption of par- 
allel incidence simplifies the analysis considerably 
and permits an analytical solution to the image re- 
construction problem; it also allows the analytical 
inversion of the complete Mossbauer spectrum at 
each point (section 3). Moreover, the parallel- 
incidence case also permits the closed-form 
derivation of the image point spread function, 
which is defined as an image of a point object 
(Dirac delta function) and characterizes the resolv- 
ing power of the imaging technique (section 4). In 
section 5 we consider an extended source which 
gives rise to y-rays incident over a range of direc- 
tions that are no longer parallel. In the latter case, 
the image reconstruction cannot be inverted ana- 
lytically as for parallel incidence, but the recon- 
struction can nonetheless be performed using 
iterative algebraic reconstruction techniques [4,5]. 
Although such iterative reconstruction methods of- 
ten fail to provide the insight and intuition of ana- 
lytical solutions (which permit, for example, the 
closed-form derivation of the image point spread 
function), such methods are often more flexible in 
incorporating a priori information into the inver- 
sion algorithm. An important example is the incor- 
poration of y-ray attenuation due to a variety of 
scattering and absorption (resonant and non- 
resonant) mechanisms. For relatively thick objects, 
attenuation can no longer be ignored in the 
inversion. For the sake of tractability, we will, 
however, ignore attenuation in the next several 
sections. 



2. Standard Experiment 

Suppose a y-ray emitted by a moving source is 
incident on a two-dimensional stationary absorber 
(fig. 2). For our purposes, it will suffice to describe 
the resonant absorption classically by modeling the 
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Figure 2-Two-dimensional absorber and distant moving source 
providing parallel y-ray illumination. 

nuclear resonance phenomenon, within a single ab- 
sorbing nucleus, as a harmonic oscillator character- 
ized by resonance frequency o>, and linewidth r 
(analogous to the mean lifetime of the y-ray energy 
level) and driven by the electromagnetic field of 
the incident y-ray. Letting A{t) denote the re- 
sponse of this elemental oscillator, we have 



A(t)+k4(t)+o>U(t)=AUt), 



where 



AjJt)=A exp /o) f 1 -^ V 



(1) 



(2) 



is the y-ray electromagnetic field. In eq (2), v, is the 
velocity of the moving source, c is the velocity of 
light and A is a constant. Although one might also 
choose to include a finite linewidth in the y-ray 
field A mc (t), we shall, for mathematical simplicity, 
lump all linewidth effects into the parameter r in eq 
(1). Note that &) v,/c in eq (2) is just the classical 
Doppler shift due to the motion of the source. In 
an equivalent notation, o) —Eo/fi and r=^/r, 
where E is the energy of the incident y-ray and T 
is the (energy) linewidth. Now in eq (1), set 



w,=o) fl-^j, 



(3) 



where v, is a parameter (in units of velocity) that 
allows for shifts in the nuclear resonance of a given 
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nucleus. Thus, the "intrinsic velocity," v„ is char- 
acteristic of different absorbing nuclei (or their en- 
vironment). We will later allow v, to vary with 
position. The solution of the differential eq (I) on 
substituting eq (2) is 






where 



<u,=w (l-Jj 



(4) 



(5) 



is the Doppler-shifted source frequency. Substitut- 
ing eqs (3) and (5) into eq (4), and assuming v, and 
v,< <c, eq (4) becomes 



A(t)= 



A 



„*».' 



2fc>o (v,— >vM/c +i'/2t ' 



Classically, the intensity of the re-emitted y-ray is 
proportional to {A | 2 , i.e., 



/(v.-v,)^M| J 



1 



(2o ;(v.-v s )Wc) 2 +(1/2t) 2 ' 
or, in an equivalent notation, 

a 2 



(6) 



« v >-H£) 2 <^ 



vjcyEi+r 2 /*' 



which is the familiar Lorentzian (or Breit-Wigner) 
distribution [1,2]. 

Because the detection of the re-emitted y-ray is 
an incoherent process, we can regard eq (6) as pro- 
portional to the probability of the y-ray re-emission 
as a function of the velocity difference v=v,— y s . 
To make eq (6) a valid probability density function 
in v, normalize eq (6) to have unit area when inte- 
grated with respect to v, i.e., define 






(7) 



where 

v P 53c/2ci> T 
is the linewidth in units of velocity. From eq (7), 



r. 



/> <r, (v)dv = l. 



One can also show that 



limP w (v)=8(v), 

r r -<0 



(8) 



where 80 is the Dirac delta function. 

Now let <r(x,y, v a ) denote the resonant absorp- 
tion cross section as a function of space and the 
intrinsic velocity v a . Then the expected number of 
counts, denoted by / (T) (v s ), is proportional to the 
function cr(x,y, v„) weighted by the probability 
density function P(v a — v s ) and integrated over v a , 
and all space: 



/ (T) (v 8 )=JJJdxdjdv a o-(x,>-,v a )/'«(v a -v s ). 



(9) 



In this and subsequent integrals, we shall assume 
that the limits of integration are from — oo to + oo 
unless otherwise indicated. We can assume this by 
defining functions such as cr{x,y,v^) to be zero 
outside their support (i.e., beyond the boundaries of 
the absorber and the domain of v a , and so on). 
Now, for convenience, suppose we temporarily 
consider the limit of very small v T (narrow resonant 
linewidth). As v T goes to zero, the probability den- 
sity function P (T) (v a — v s ) becomes highly peaked 
about v a =v s , and in the limit v T —*0, we can use eq 
(8) in (9). Thus, defining 



/(v^lim/^v,), 



we have from eqs (8) and (9), 

/(v s )= I I dxdyor(x,y, vj, 

which is the expected result for the standard Moss- 
bauer experiment, that is, /(v s ) is the absorption 
cross section evaluated at v a =v s averaged over the 
sample. 

3. Imaging with Parallel y-Ray Incidence 

Now let the absorber rotate clockwise about the 
z-axis at constant angular velocity ft. Let (x,y) de- 
note a stationary coordinate system and let 
(x',y') denote a moving coordinate system embed- 
ded in the rotating absorber, as illustrated in figure 
3. At r=0, let the stationary (x,y) and the rotating 
(x',y') coordinate systems coincide; at this instant 
denote the absorber's cross-section function by 
°"o(*t.)', v a ). For r>0, the absorber has rotated 
through an angle £lt. Define cr,(x,y, v a ) as the ab- 
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sorption cross section relative to the stationary axes 
(x,y) at time t; then 



<r,(x,y, v a )=or (x',^', v a ), 

where 
x'=x cos ilt — y sin fit 
y'=x sin fit +y cos ftf , 

or, inverting, 
x — x ' cos fir +y ' sin ftf 
y=—x' sinClt+y' cos fir. 



(10) 

(11a) 
(lib) 

(12a) 
(12b) 



Again consider a y-ray incident parallel to the 
x-axis from a source moving at velocity v s in the 
positive x -direction (fig. 3). Resonant absorption 




-fh 



Condition lor resonant absorption: RQ = v s -v a 



Source (distant) 



Figure 3-Rotating two-dimensional absorber with parallel y-ray 
illumination. Resonant absorption takes place only along a 
line of "sensitivity" parallel to the x-axis at height R, where 
R is defined by the condition Ril = v t — v,. 

will then take place when the total velocity v,+v,' 
of a point on the absorber equals the source veloc- 
ity v s , where v a is the intrinsic material "velocity" 
and v„' is the additional velocity contribution due to 
the motion of the absorber. Since the y-ray is inci- 
dent parallel to the x-axis, v a ' will be the x -compo- 
nent of the absorber velocity; thus, from eq (12), 

v a '=x=fij>. 

Note that v a ' depends only on y, so that resonant 
absorption takes place along lines parallel to the 
x-axis. For the stationary absorber, the resonance 
condition was v a =v s ; for the rotating absorber, the 
resonance condition is modified to read 



Thus, for specific values of ft and v„ resonant 
absorption takes place along a line at height 
y=(v s — v,)/ft parallel to the x-axis (fig. 3). 

Now, in view of eq (13), for the rotating ab- 
sorber, we can replace v, in eq (9) by v,— fly, giving 



/< r >(v s , ft)= (((dxdydv t o;(x,y, v.) 
.pV(v,-v.+Sly). 



(14) 



Now considering a very narrow resonance, we 
again take the limit v T — >0. Defining 



/,(v„ ft)=lim/, w (v„ ft), 

v T ->0 



and using, from eq (8), 



lim P (T) (v a - v,-r-fty)=8(v.- v.+fty) 



v t +ily=v i . 



(13) 



=^S[v-(v.-v.)/ft] 



in eq (14), we have 

f,(v it il)~ JJdxdv.cr,[x, (v,-v.)/ft, v.]. (15) 

This is our fundamental equation from which we 
shall solve for the function crd(x',y', v.) from the 
data/,C.» ft) recorded as a function of v„ ft, and /. 
Below, we consider four special cases of increasing 
generality: 1) one-dimensional imaging with no v, 
distribution (i.e., assuming a single spectral line at 
v,=0); 2) two-dimensional imaging with no v, dis- 
tribution; 3) one-dimensional image with an arbi- 
trary and unknown v, distribution; and finally 
4) two-dimensional imaging with an arbitrary and 
unknown v, distribution. 

One-Dimensional Imaging with No v t Distribu- 
tion. We consider the simplest imaging case here. 
Suppose the absorber is one-dimensional and rotat- 
ing. For convenience, assume the absorber is 
aligned with the (rotating) >> '-axis, and suppose fur- 
ther that the measurements arc made at the instant 
the absorber rotates past vertical, i.e., when the ro- 
tating absorber coincides with the y-&xis at f =0. 
We then have 

tr,{x,y, v,)=cro(y)8(x)6(v,). 

and eq (15) reduces to the simple result 
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Thus, the spatial dependence of the one-dimen- 
sional cross section cr (y) is given explicitly as a 
function of vj(l. 

Two-Dimensional Imaging (Tomography) with No v a 
Distribution. Here we let 



<r,(x,y, vj=cr t (x,y)8(v,), 
so that eq (15) becomes 

f,(v„n)~jdx<r,(x,v l /il). 



(16) 



(17) 



Recall that the measurements/^,, ft) are made as a 
function of time as the absorber rotates. To make 
the /-dependence in eq (17) explicit, rewrite eq (17) 
as follows. 

/(»„ ft)=^ j(dxdycr,(x,y)8(y-v i /(l) 

=^jy&c'd/0r (x',/) 

•S(-x'sinnr+^'cosn/-v/n), (18) 

where eqs (10) and (12b) were used in the last step. 
From eq (18),/(v„ ft) is a line integral through the 
two-dimensional function <r (x',y') along a line 
parameterized by the radial distance v/ft and the 
angle (It from the .y'-axis. The inversion of eq (18) 
for <rd(x',y') can be carried out by standard tomo- 
graphic techniques. It is easiest to invert eq (18) 
using Fourier methods, as follows. 

First define the one-dimensional Fourier trans- 
form of/,(v„ H) with respect to v, as 

/,(*, ft)= jdv,/(v„ ft)exp(-/Xv 5 ). (19) 

Inserting eq (18) into (19) and interchanging orders 
of integration, gives 



■/;■ 



/,(*.ft)=JJdx'd/cr (jc',j.') 

•exp[i'&n(*' sin (It -y' cos (It)). (20) 

Now the two-dimensional Fourier transform, 
d-<,(k„k r ), of a (x ',y') is 



o-o(*„ *,)= IT dx'dy'<ro(x',y'. 



) 



•exp[-/(^'+V)], (21) 

and, on comparing eqs (21) and (20), we have 

f t (k, ft)=6- (-/fcft sin fir, k(l cos (It). (22) 

Equation (22) is the well-known central-slice 
theorem in tomography [4,5]. In the present con- 
text it states that f,(k, ft), evaluated at time t and 
spatial frequency k, is the two-dimensional spatial 
Fourier transform, & (k x , k y ), of cr (x',y') evalu- 
ated on a line in the (k x , k y ) Fourier plane through 
the origin at angle (It with the k x -axis. Thus, taking 
many revolutions of data and, during each revolu- 
tion, stepping the value of v s will generate sufficient 
data to provide complete coverage of the 
(k x , k y ) Fourier plane, whereupon an inverse two- 
dimensional Fourier transform of cr (&,, k y ) yields 
the reconstruction cr Q (x',y'). Note that the re- 
quired range of measured v s values is from to 
ft^maw where /? mM is the largest radial dimension of 
the absorber. 

One-Dimensional Imaging with Unknown v a Distri- 
bution. In this case, let 

<r,(x,y, Vi)=o- (y, v*)8(x). 

For the one-dimensional absorber, again assume 
that the measurement is made as the absorber ro- 
tates past vertical at r=0. Then, substituting 
<r,(x,y, y a ) into eq (15) and setting t=0, gives 

Mv u ft)=^ JdVaO-oKn-vJ/ft, vj. (23) 

We wish to solve for the two-dimensional function 
o-oOs vj, which has one spatial dimension (y) and 
one spectral dimension (vj. Rewrite eq (23) as 

f&„ ft)=^ JTdj>dv,o-oO'. v a )8[y -(v.- v a )/ft]. 

Substituting into eq (19) and interchanging orders 
of integration, gives 



/, 



(k, ft) = d^dv a cr o 0;, v a )exp[-/£(ftj> + v a )] 
J J i 



(24) 



Now writing the two-dimensional Fourier trans- 
form of o- 0\ y .)> 
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ar {k y , fc v )= dydv a o- ()>, K 



) 



•exp[— i(Jc,y+k r v^], 

and comparing to eq (24), we see that 

fo(k,ii)=&o(.ka,k). 

This result resembles the central-slice theorem en- 
countered earlier. Here f (k, ft) is equal to the 
two-dimensional Fourier transform, & (k y , k,), of 
o r o(y» v») evaluated on a line in the (k„ k t ) Fourier 
plane through the origin at angle tan~'(ft) from the 
fc,-axis. Thus, by varying both k and ft, complete 
two-dimensional coverage of & (k y , k„) in Fourier 
space can be achieved. The reconstruction of 
o"o(y, v a ) then follows on taking the inverse two-di- 
mensional Fourier transform of & Q (ky, k,). 

Two-Dimensional Imaging (Tomography) with Un- 
known v, Distribution. Here we wish to recon- 
struct the three-dimensional function cr (x',y',vj. 
This is the most general inversion problem consid- 
ered. Now rewrite eq (15) as 

/(v s , ft)=jj dxdydv a cr,(x,y, v a ) 

•8[y-(v s -v a )/ft] 

= ft 1 1 1 to'd/dv.oroCx',/, O 

•8[-x' sin Sit +y' cos ftf -(v s - v a )/ft] 

(25) 

using eqs (10) and (12b) in the last step. Substitut- 
ing eq (25) into eq (19) then gives 

f,(k, n)= f f fdVd/dv a o- (x',/, v a ) 

•exp[/A:ft(x' sin ftf —y' cos ftf)— /JfcvJ. 

(26) 

Writing the three-dimensional Fourier transform of 
°o(x',/, v a ) as 

o-o(.k x , k y , fc v )=J (J dx'dy'dvtcr a (x\ v', v.) 

•exp[-f(M'+fc^'+fc,v.)], 
and comparing to eq (26), we have 



f,(k, ft)=er (-fcft sin fit, k£l cos ftf, k). 

This is the generalization of the central-slice 
theorem to three-dimensions. Thus, by suitably 
varying the three parameters v„ ft and /, full 
coverage in the three-dimensional Fourier space 
can be achieved, and an inverse Fourier trans- 
form of & {k x , ky, k,) yields the reconstruction 
cr^x',y', v.). 



4. Spatial Resolution and the Image Point 
Spread Function 

In previous sections, we considered the idealized 
case of an infinitesimally narrow linewidth v r . This 
simplifies the mathematics, but implies that spatial 
resolution can be arbitrarily small. In this section, 
we show that spatial resolution is controlled by the 
natural linewidth v T ; in particular, spatial resolulion 
is proportional to the ratio of v T to the rotational 
velocity ft of the absorber. 

Below we derive the point spread function 
(PSF) of the two-dimensional tomographic recon- 
struction problem, with no v, dependence for sim- 
plicity (the second case above). The PSF is by 
definition merely the reconstructed image of a two- 
dimensional Dirac delta function, or point object, 
and its width provides a reasonable measure of im- 
age resolution. Alternatively, the reconstructed im- 
age may be regarded as the true image convolved 
with the PSF. 

To derive the point spread function in two di- 
mensions, we repeat the above derivation using eq 
(14) in place of eq (15). Substituting 

(r,(x.V.V.) = 0r,(x,v)8(v.) 

into eq (14) gives 

fl>\v t ,n)=jjdxdyo;(x,y)P«\ny-v.), 

which is the first line of eq (18) with §(•) replaced 
by P <r) (0- Now write this as 



^{ V „il)=jjdx , dy'o-^x',y') 



•P ( "[ft(-x' sin ftf +/ cos ftf)- v,]. 

(27) 

Writing the Fourier transform of/T'fv,, ft) with 
respect to v„ 
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/«(*, a)=(dv t /}*(r„ ft)exp(-/*v s ), (28) 

substituting eq (27) into (28) and interchanging or- 
ders of integration, results in 

/«(*, n)=P"(k)jjdx'dy'cro(x',y') 

•exp[/&ft(*' sin ft/ —y' cos ft/)] 

= P (T, (A:)cro(-A:ft sin ft/, &ft cos ft/), 

(29) 



where 



P (r, (*)= [dv s P w (- v s )exp(-/A:v s ) 



(30) 



is the Fourier transform of P (r, (— v,). Equation (29) 
is merely eq (22) multiplied by P (T) (fc). 

We can evaluate P in (k) by inserting eq (7) into 
(30); integrating, the result is 

P <T \k)=exp(-kv T ). 

Thus, the Fourier transform, & (Kt k y ), of 
cr (x',y') is weighted by the exponential radially- 
dependent function exp(— kv r ). The point spread 
function is obtained by setting &d(k x , k y )=l in eq 
(29); as a result, the PSF is just the two-dimen- 
sional inverse Fourier transform of P {T) (k). Since 
kCl is the radial coordinate in the Fourier plane, 
P (T) (k) has radial symmetry and the resulting PSF 
also has radial symmetry. Hence, the inverse two- 
dimensional Fourier transform of P ( '\k) can be 
evaluated in polar form: 



ft 2 f" 
PSF(0=|J fl 



dk k P< T) (k)Jti(nkr') 



=~ f dkk e\p(-kv r )J ((lkr') 



Vr 



1 



(31) 



^KO'+^/ft) 2 } 3 ' 1 * 

where J (-) is the zero-order Bessel function and 
r'=\/x' 2 +y' 2 . This is the desired point spread 
function whose width gives the image resolution. 
Let r ' signify the full width at half maximum of eq 
(31); then 

r '=0.76v r /ft. 

The important thing to note is that r ' is inversely 



proportional to the rotational velocity ft. Finally, 
one can show that PSF(r') becomes a two-dimen- 
sional Dirac delta function in the limit v T ->0, as 
required. As an example, for typical linewidths, v,, 
of less than a mm/s, this result shows that r ' is on 
the order of a tenth of a millimeter when the ab- 
sorber rotates at one revolution per second 

(ft = 277). 



5. Imaging with an Extended Source 
(Non-Parallel Incidence) 

In the latter two sections, we assumed parallel 
y-ray incidence in the direction of the a: -axis. Par- 
allel incidence is easily achieved by placing the 
source far enough away so that it subtends a small 
angle at all points on the absorber. This is practical 
when the source is sufficiently strong to provide a 
■y-ray flux intercepting the absorber that results in 
an acceptable signal-to-noise ratio. Unfortunately, 
arbitrarily strong Mossbauer sources are not avail- 
able, so that moving the source close to the ab- 
sorber to maximize the incident flux may be 
important. 

In the latter case, the source must be regarded as 
spatially extended, say, in the j> -direction (fig. 4), 
and the incident y-rays are no longer parallel. As- 
sume that the extended source is moving in the x- 
direction at velocity v s and that the absorber rotates 
with angular velocity ft as before. Also, for sim- 
plicity, assume that the intrinsic velocity v a =0. In 
this geometry, one can show that the component of 
velocity of a point on the rotating absorber and of 
a point on the source in the direction of the line 
joining them (i.e., along the line of flight of a y-ray 
passing between the two points) are equal to each 
other (the resonance condition) if and only if the 
line defined by the two points passes through 
(x=0, y = vj£l) in the stationary (x,y) coordinate 
system. To see this, refer to the line L in figure 5, 
which we suppose is the path of a y-ray emitted 
from point A on the source. Suppose this path in- 
tersects the j'-axis at height R , as shown. Then L is 
tangent to a circle around the axis of rotation of 
radius R cos 0, where is the angle between L and 
the ;r-axis. From previous arguments, all points in 
the absorber lying on L have the same component 
of velocity in the direction of L, namely fli? cos 6. 
But the additional component of velocity of the 
incident y-ray along L due to the motion of the 
source is v s cos 0, as can be seen from figure 5. The 
resonance condition is then obtained by equating 
the absorber and source components of velocity 
along L; thus setting HR cos 0=v s cos results in 
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R =Vs/fl. Consequently, only those y-rays moving 
along paths through (x =0, y = v/fl) are resonantly 
absorbed. 

The complete set of lines passing through this 
point and intercepting the source generates a sec- 
tor, as shown in figure 4, where the boundaries of 
the sector are the two lines intersecting the end 
points of the extended source and coming together 
at the point (0, v/ft) in the (x,y) system. The den- 
sity of lines peaks at the point (0, v/fl) and falls off 
in proportion to l/\x\ in the x -direction (fig. 6). 
Through the peak in the ^-direction, the width of 
this function is much narrower; in particular, the 
width in the y -direction can be shown to be v/ft, 
which is approximately the width of the point 
spread function derived previously for parallel in- 
cidence. Neglecting, for the moment, attenuation 
of the y-rays through the absorber, the line-density 
function sketched in figure 6 may be equivalently 
interpreted as the probability density function of 
resonant absorption events per unit area, which is 
just proportional to the inverse of the sector area 
shown in figure 4. 

The highly peaked distribution of the probability 
of resonance events is significant, since it implies 
that merely by varying v s and counting absorption 
events as a function of time, it should be possible to 
generate a blurred image even without image re- 
construction. Figure 4 shows that, by rotating the 
absorber, the y-rays are effectively "focused" 
through the point (0, vJCl). This focusing-like ef- 
fect can be used to create a blurred image without 
further processing. A complete or deblurred image 
reconstruction cannot in this case be carried out 
analytically, but the reconstruction nevertheless 
can be performed using iterative algebraic recon- 
struction techniques (ART) [4,5]. The ART al- 
gorithms assume an initial estimate of the unknown 
image and then, on the basis of this estimate, re- 
compute the measurements (line integrals of the 
resonant absorption coefficient). The recomputed 
measurements are then subtracted from the true 
measurements, giving rise to a set of error terms, 
which are then used to systematically adjust the 
current image estimate to reduce this error. When 
this iterative procedure is continued, it can be 
shown that the overall error tends to zero and, in 
the absence of measurement uncertainty and in- 
complete data, the reconstructed image converges 
to the true image. Such iterative approaches also 
make it relatively easy to incorporate attenuation 
effects in the algorithm. 

In the above discussion, we assumed for simplic- 
ity that the extended source was one-dimensional 
(i.e., extended along the y -direction). If the source 



has some depth as well as height (i.e., some exten- 
sion in the z-direction), the additional degree of 
freedom arising from variations in the x -compo- 
nent of the y-ray velocity will broaden somewhat 



Sector of 
resonant 
absorption 




Figure 4-For an extended source (non-parallel y-ray incidence), 
the line of sensitivity broadens into a sector of sensitivity 
bounded by lines subtending the source and intersecting at the 
point (x=0,y=R), where Rtl=v,-v,. 
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Figure 5-The line L signifies the path of a y-ray emitted from 
one point on the extended source. 




Figure 6-Probability distribution of resonant absorption per unit 
area over the sector of sensitivity in figure 4. The distribution 
falls off as l/|jt| along the x-axis and is approximately of 
width v r /fl along the .y-axis. 
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the region of sensitivity beyond the sector illus- 
trated in figure 4. To minimize this broadening ef- 
fect, the source should be collimated in the 
z -direction as much as practical with a suitable ar- 
rangement of slits in this direction. A three-dimen- 
sional absorber could then be reconstructed on a 
"slice-by-slice" basis, as in two-dimensional x-ray 
tomography. 



6. Applications 

Applications of the technique to interior imaging 
appear to be limited to lighter materials, or ab- 
sorbers of small size, that permit sufficient y-ray 
penetration. One possibility is the imaging of com- 
posite materials, in which a Mossbauer element is 
embedded in a lighter matrix. Another possibility is 
surface imaging, whereby the surface of a flat ab- 
sorber is rotated about an axis normal to the surface 
and a one-dimensional extended source is displaced 
a small amount from the plane of the surface. If all 
the impinging -y-rays have the same (or nearly the 
same) z -component of velocity, the region of sensi- 
tivity in which resonant absorption takes place will 
look much the same as the sector illustrated in fig- 
ure 4. In this case, the value of v, used previously 
should be replaced by the component of the source 
velocity in the plane of the surface. 

Some potential high-resolution applications in 
materials science include the imaging of grain 
boundary segregation, the imaging of residual 
stress distributions, and the imaging of the distribu- 
tion of magnetism in ferromagnetic materials due 
to perturbations in the Mossbauer spectra arising 
from variations in the magnetic hyperfine field. 



image of the absorption coefficient distribution in 
the absence of further processing when the mea- 
surements are recorded as a function of v s and time 
(or the instantaneous absorber rotational position). 
Furthermore, by exploiting the three measurement 
degrees of freedom, v s , ft, and time (or rotational 
position), complete spectroscopic information can, 
in principle, be recovered as a function of position. 
Thus, Mossbauer spectroscopy can be performed 
in an imaging mode. Spatial resolution is propor- 
tional to the ratio of the natural linewidth v T to the 
rotational velocity 11, but becomes rapidly signal- 
to-noise limited as the resolution increases. 
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7. Conclusion 

In conventional Mossbauer spectroscopy, the 
measurement is spatially integrated over the ab- 
sorber. For parallel y-ray incidence, a rotating ab- 
sorber generates a line of "sensitivity" over which 
resonant absorption takes place, where the location 
of the line is a function of the source velocity and 
the absorber rotational velocity. In particular, this 
line is parallel to the velocity of the source and 
passes through the point (x=0, y = v/fl). For an 
extended source, the line broadens into a sector 
subtending the source and coming together at the 
point (x=0, y = v/Cl). For a highly extended 
source, the result is as if the y-rays were effectively 
"focused" through the point (x=0, >>=v,/ft). This 
focusing effect can be used to generate a blurred 
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